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A multiple-moment approach to the dielectric function of a dense non-ideal plasma is treated 
beyond RPA including collisions in Born approximation. The results are compared with the pertur- 
bation expansion of the Kubo formula. Sum rules as well as Ward identities are considered. The 
relations to optical properties as well as to the dc electrical conductivity are pointed out. 



I. INTRODUCTION 

The dielectric function e(fc, lo) is a physical quantity containing a lot of information about the plasma. In homo- 
geneous, isotropic systems it is related to the electrical conductivity a{k^uj) and the polarisation function T\{k,uj) 
according to 

e(fc, uo) = l + — <j{k, uj)^l ^ n(fc, uj). (1) 

A well established expression is the random phase approximation (RPA) valid for coUisionless plasmas. The inclusion 
of collisions, however, is connected with difhculties. A perturbative treatment of the Kubo formula is not applicable 
near fc = 0, w = because there is an essential singularity in zeroth order. Partial summations are sometimes 
in conflict with sum rules. Improvements of the RPA result are discussed in the static limit, where the local field 
corrections are treated in time-dependent mean-field theory Also approximations based on the sum rules for the 
lowest moments have been proposed Q. However, an unambigous expression for e(fc, lu) in the entire k lo space cannot 
be given by these approaches. 

A particular problem is the appropriate treatment of the long-wavelength limit fc — > at small frequencies where 
the dc conductivity should be obtained. In a previous paper an approach has been given where this limiting case 
coincides with the Chapman-Enskog approach Q to the dc conductivity. In particular, the polarization function was 
found as 

n(fc,tj) = i — /3rio 

LO 

The matrix elements Mmn are equilibrium correlation function which are explicitly given in the following section. 

They contain operators Bm and i3„ which specify the nonequilibrium state. 

For the evaluation of the dielectric function, we have to deal with two problems: 

i) the choice of the operators to describe the relevant fluctuations in the linear response regime, 

it) the evaluation of the equilibrium correlation functions. 

The equilibrium correlation functions in a nonideal plasma can be evaluated using the method of thermodynamic 
Green functions. In lowest order of the perturbation theory to be considered here we have the Born approximation 
as described in . Higher order terms can be taken into account in a systematic way, see Q . 

With respect to the choice of the operators i?„ , only the current density operator J has been considered in Q . In the 
spirit of the Chapman-Enskog approach we will include here higher moments of the single-particle distribution function 
to study the convergency behavior. For the dc conductivity the answer is well known see 0. Note that different 
approaches based on different sets of relevant observables i?„ are formally equivalent as long as no approximations in 
evaluating the correlation functions are performed. However, within a finite order perturbation theory, the results for 
the conductivity are improved if the set of relevant observables is extended. 

Results for the dielectric function within a four-moment approach are shown in Sec. II and compared with the 
results of a single-moment approach. Some exact relations are discussed in Sec. HI. Of particular interest is the 
relation to the Kubo formula which may be treated in perturbation theory as discussed in Sec. IV. 



Mo„(fc,t^) 

M^Q(k,Lo) Mmnik,(o) 



\Mrnn{k,L0)\. 



(2) 
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II. FOUR-MOMENT APPROACH FOR A TWO-COMPONENT PLASMA 



To evaluate the dielectric function we use the expression (H) for the polarisation function, where the matrix elements 
are given by 



Mmn(k, LO) = [Bm] [Bn - i^Bn]) + [S„ - iuj B^]) uj+ir, - ,„ ' , {Bm] [-B„ - iuj Bn]) uj+ir, ■ (3) 

The equilibrium correlation functions are defined as 



{A; B) = {B+;A+) ^\fdT Tr [A{~inT)B+ p^] 
P Jo 

{A;B)^= r 
Jo 



dte^^'{A{t);B) , (4) 

withyl(<) = cxp{iHt/h) A eyii>{-iHt/h) and i = . po = eiip{-(3H + f3J2c f^cNc) j Trexp(-/3i/ + ^^^/icA^c) 

is the equilibrium statistical operator. 

We will consider a two-component plasma consisting of electrons (c — e) and ions (c = i). In particular, results are 
given below for a hydrogen plasma. With the single-particle operators 



n. 



the current density operator is given by 



c.p 

Furthermore we used the abbreviation = (r^[k,uj) Jk . 

To select the relevant operators Bn, we restrict us to the ordinary kinetic approach. The inclusion of higher order 
correlations is also possible, see 

Within the kinetic approach, the nonequilibrium state of the plasma is described by the mean values of th single- 
particle operators corresponding to an induced single-particle distribution function with wave number k. Instead 
of treating an infinit number of operators depending on the momentum p, we can restrict us to a finite number of 
moments of the distribution function. This procedure is familiar from the theory of the dc conductivity. Whereas in 
that case only moments with respect to p have to be selected, in the general case of arbitrary k to be considered here 
moments of p as well as p- k have to be taken into account. 

In this paper we investigate how the lowest moment approach in Born approximation is modified if further moments 
are included. From the theory of dc conductivity we know that important modifications are obtained by including 
the energy current density in addition to the particle current density, i. e. if we include also "j^Pz- Then, the electrical 
cunductivity is not only described by the electron-ion interaction, but includes also the effects of electron-electron 
interaction which are not effective in the lowest moment approximation due to the conservation of total momentum. 

The four-moment approach to be considered in this paper is given by the following moments of the electron (c = e) 
or ion (c — i) distribution function, respectively, 

\/ Zmckl 

The evaluation of the corresponding correlation functions in Born approximation is given in the Appendix for the 
nondegenerate case. As a trivial result, in the lowest approximation with respect to the interaction the RPA result 
is recovered. In general the matrix elements are given in terms of integrals of expressions containing the Dawson 
integral. 

To give an example, a hydrogen plasma is considered with parameter values T = 98 Ryd and rie — 8.9 a^'^ which 
are found in the center of the sun |7]. The results are comparable to the results obtained in for parameter values 
corresponding to laser produced high-density plasmas S. 
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Results for the real and the imaginary part of the dielectric function in the two-moment approximation given by 
bKp) are shown in figures 1 and 2, respectively. Besides the RPA dielectric function the one-moment calculation 
reported in [|j is shown as well. While the differences between the improved dielectric function and the RPA are small 
at high momenta (k = 1 a^^ ), significant changes occur at small momenta (fc = 0.1 cig^)- On the other hand, the 
one-moment approach is almost identical with the two-moment calculation. This is an indication that convergence is 
reached by augmenting the number of moments as is expected from earlier studies of the dc conductivity Note, 
that the static limit is given by the Debye law. 

Results for the inverse dielectric function, which describes the response to the external potential, are shown in 
figures 3 (fc = 0.5 a^^) and 4 {k = 0.3 a^^) and compared with the RPA inverse dielectric function. Major deviations 
occur only at frequencies close to the plasma frequency. For small momenta, the imaginary part of the dielectric 
function including collisions is considerably broader compared with the RPA one. While the imaginary part of the 
inverse dielectric function in the RPA approximation becomes delta-like in the long wavelength limit, a broadening of 
the plasmon peak appears, as can be seen from figure 4. Some properties of the dielectric function will be discussed 
in the following section. 



III. EXACT RELATIONS FOR THE DIELECTRIC FUNCTION AND LIMITING CASES 

Several exact properties of the dielectric function are known such as sum rules 

duj 



the long-wavelength limit 



ujlme ^(fc,w) = ujpi , (8) 

oo TT 

— wlme(fc,cj) = , (9) 

-OO 

/OC 1 1 
— -lme-\k,iu) = -l , (10) 
-oo I" ^ 

limRee(fc,0) = 1 + F(fc)n2if . (11) 



fe^O 



Here ~ X]c=e i (^^ ^c) / (eo ^c) denotes the plasma frequency and K the isothermal compressibility. Further 
extensions for a two-component system can be found in ||l^. This is a special relation resulting from the relation 
between the dynamical structure factor 

S{k,Lu) = — / dt <p+{t)pk> e^"^* (12) 



2 7r „ 

and the dielectric function which can be established via the fluctuation-dissipation-theorem 

S{k,oj)^--^^lme-\k,oj-) . (13) 

Furthermore, the Kramers-Kronig relation holds which connects the real and the imaginary part of the dielectric 
function: 

Ree{k,uj) = 1 + P ■ 14 

J n u) — u)' 

Here, P denotes the Cauchy principal value integration. The inverse dielectric function obeys a corresponding relation. 
Combining the Kramers-Kronig relation with the sum rules results in rigorous statements about the asymptotic 
behaviour at high frequencies: 

lim Ree(fc,cj)==l-^^ + 0(^) . (15) 

UJ— *oo iO UJ 

We test the two-moment approach by checking the sum rules as well as the asymptotic behaviour. It is found 
that the sum rules are fulfilled within the numerical accuracy (« 0.1%). The Kramers-Kronig relation holds as well. 
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Having in mind relation the dc conductivity can be obtained considering the hmitng case k ^ 0. A comparison 
with other results for the conductivity can be made by parameterising the conductivity via 

a = a(0,0) = s — , (16) 

where $ denotes the Coulomb-logarithm and rriei the reduced mass of the electron. 

As shown in figure 5, there is no shift of the maximum of the inverse dielectric function, while the plasmon peak is 
broadened. Moreover, the long wavelength limit can be described by a Drude-like formula, implying that the width 
of the plasmon peak is given by the dc conductivity. The form of the plasmon peak can be compared with computer 
simulation studies. In contrast to RPA calculations, width as well as height of the plasmon peak in our calculation 
are of comparable size as computer simulations [p^ . 



IV. COMPARISON WITH THE KUBO FORMULA 



Depending on the selected set of relevant operators {Bn}, different expressions for the dielectric function can be 
derived within linear response theory. A often used expression is the Kubo formula jist as given by 

n(fc,w) =- {Jk.Jkluj+ir, ■ (17) 

UJ 

As shown in this result follows as a special case within the generalized linear response theory. As also shown 
there, the different expressions identical in the limit 7y — > if no further approximations are performed. 

The advantage of linear response theory is that the evaluation of the dielectric function is related to the evaluation 
of equilibrium correlation functions. In dense, strongly coupled systems, these correlation functions can be calculated 
with computer simulations. Another possibility is to use peturbation theory which is most effectively formulated with 
the concept of thermodynamic Green functions [ p^ . 

In zeroth order with respect to the interaction, from ( [l7| ) immediately the RPA result is obtained, in coincidence 
with all other approaches including J within the set of relevant operators. The first order expansion with respect to 
the screened interaction reads 



n(fc, ^,) = E ( /p + E -E ) 



1 1 



Ep- uj\- Ep-k Ep + ujx- Ep+k , 
kq 1 1 1/1 1 



m Ep — Ep^q Ep — uj\ — Ep^k Ep — ux — Ep^k-q \Ep — Ep^g Ep — uj\ — Ep^k , 
+ (a;,fc ^ -w,-fc). (18) 

For the sake of simplicity, we have taken the adiabatic limit where rrii/mf, — > oo (Lorentz plasma). In particular we 
find for fc ^ 

lmn{k,u) = nY,V^(^\ 7rS{Ep-uj-Ep^g)e'P^^--^^^-^^. (19) 



m 

pq 



what gives the frequency-dependent conductivity. 

However, this perturbation expansion does not converge at — > 0, and partial summations have to be performed. 
For instance, a simple approximation for the polarization function including interactions with further particles would 
be a polarization function given by the product of two full propagators. This way, the polarization function contains 
shifts and damping of the single-particle states due to the interaction with the medium. However, this approximation 
does not fulfill rigorous relations such as sum rules, since important corrections to the RPA of the same order in the 
density as the considered ones are missing, e.g. vertex corrections. These corrections are linked to the self-energy by 
Ward identities ||l^. As a consequence, the vertex has to be improved in accordance with the self-energy. Following 
Baym and Kadanoff ||l^, a consistent vertex can be constructed to a given self-energy. However, the solution of the 
vertex equation cannot be given in a simple algebraic form, and usually some approximations are performed, see [ pT| . 
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V. CONCLUSIONS 



An approach to the dielectric function has been investigated which includes the effects of coUissions and can be used 
in the entire k,uj space. Within a four-monicnt approach to a two-component plasma, the Born approximation has 
been evaluated, and important rigorous results for the dielectric function are checked. Compared with the ordinary 
Kubo formula, the approach given here seems to be more appropriate for perturbation expansions. 

In particular, comparing with a one-moment approach, the convergency behavior of this method was inspected. As 
well known from the theory of dc conductivity, convergence is expected if higher moments are included. In a more 
general approch, also two-particle correlations can be included into the set of relevant operators. 

Within a quantum statistical approach, the Born approximation can be improved by systematic treatment of Green 
functions. This concerns, e.g., the inclusion of strong collision by treating T- matrices, degeneracy effects, and the 
treatment of the dynamic screening of the interaction. Here, the comparison with computer simulations is also an 
interesting perspective. Work in this direction is in progress. 



The authors acknowledge helpful discussions with W. Ebeling, W.D. Kraeft, D. Kremp, R. Redmer, and Chr. 
Toepffer. 

APPENDIX: EVALUATION OF THE MATRIX ELEMENTS OF n 

We start from the general expression (1), (2) for the dielectric function with (n, m =1...4, c, = e, i) 

MLM = 7^ E ^-,^PMk){ni-,,n;,^) , (20) 

MUd,^) = 7^ E ^<i^^UK{p)rK-,^ni^,) , (21) 

a,p,k 



A^;^i(9,^) = 7j-E[^n(^')^^^^^(^)(K-«-^'^<-,];«;,,)+^ (22) 

or, after some rearrangements, 

""'^ ° -mtMi (23) 



e{q,(j) = 1 



eoqui 



with 



MLM = ^M^mM = ^7^ E ^^^•P^&™(fc)«-,;";,,) ' (24) 



c,p,k 



d,p,k 



Mt^iq^u:) = -i^^^Mi,i,{q,u) = -^^^^l-J2K{p)rht{k) {([<_, - *^<-J;n^,,) 



, l-d r-c ■ c l\ ("-fc.-qi "^fe)a>-FiTi , ^ -en 1 lc-,r\ 

+ K,g - I'^np^qUuj+iv - I d T\ ("/c,-«; K,<Z - ) , (26) 



and 
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We specify to a four-moment approach where Bi = bl{p), B2 = ^'2(^)1 -^3 — ^i(-P)' -^4 ~ ^i(-P)- Introducing 
the Dawson integral 

= lini / da:e--' (28) 



and using the abbreviations 



c 1 



1 2 l + z,D{z,) 



we have for = Ui = n, Ce = —Si = e 



(29) 
(30) 



Ki^l, K3 = -l^ ^04 = -J, (31) 



and 



We decompose 



Mfo-^, M!o = l, ^^30 = -^, Mlo--^, (32) 



M^iiiq, Uj) = a„„ + bnm + Cnm (33) 

and find in zeroth order with respect to the interaction 

le9 ^ eQ 5e9 35 

On = xri-, ai2 = -ri-, aai = -r2-, a22 = -5-7-2-, (34) 

Z Ul 4W 40^ oW 

and 

1*9 5^9 5 , g 35 , g 

033 = ;^ 7-1-, a34 = -ri-, 043 = -rj-, 044 = —r2-. 35) 

The bnm contain the electron-ion interaction in first Born approximation and the Cnm the electron-electron or 
ion-ion interaction, respectively. We use a screened interaction with the Debye screening factor exp(— Kr), = 
j:^n,el/{eokT). 

Terms due to electron-ion interaction are with M = me + m.i 

i 1 e"* / 1 \ /TTiemi \ 1/2 
- 8(2^-"-'-^ ' 

With 



and 



Amemikl p 



Ai = [In I — + — ] , (39) 
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' A"^* — 1 \ 
A2 = [A-ln(^-^ +2], 



As 



(A«^)2 - 1 



(40) 



— 



M 



(41) 



/ 771 " / 777 

De = -D(2:ei - , —pc), Di = D{Zei - \ —pc) , 

y nie y nii 



(42) 



we find 



511=^ dpe-p" Ai l^^p - Rl J rfcc(2£)e + A)| , 



(43) 



513 = ^ dpe-p" A^i^-'^p + Rl J ^dccD^j , 



(44) 



912 



I dpe-^' <jAi(-^p+-^p=^ + i?^2^ 



^3 3 



-P)+ 



+ 



J dcRl 



. , 5 TOe n _ me , 



^-pc) ) + A22p ^ ^ (1 - 3c )(^ei - 



+ AAic(--)) 



-pc) 



(45) 



ffl4 



/•oo r 



5 mi 2 me 3 ,/mJn~mi 
^^(-3M^-3M^ -"^^^ Mm. 



dcRl{De 



pc ) + A22p-!^— (1 - 3c ){Ze 

me M 





pc) 



(46) 



921 = / dpe 



5 me 2 mi 2 , oe/^ym^ mey^mJrTi 
^^^(SM +3M^ +^^^^3^1^ + ' m,M » 



dcR^{De 



+ Di 



Aic(-2 
Aic(- 



me 
M 

TOe 



TO, 2 



y/meTrii hTH me. 



'%c)^)+A22p^^(l-3c2)(^e.- 
mp iW 



3^(7 



/TOe a o y/'^eirii 2x/ 

-pc)^) - A22p-^^— — (1 - 3c^)(2;ei - 
mi Az 



— pc) 



-pc) 



(47) 



523 



dpe P Aip - 



rfciJjAAic 



5 me 
317 



m„ mi 



2 TOi 2 ne 2 a/to^TO^ 

3M^ +^23^^^ 



M M 



VTOeTO, 

M 



■P(2ei - \ —pc) + -Tj{Ze 

" me -M 



TOe 



(48) 



/■°°^ _„2 47 m2 
,22 = dpe^ f^(yMl^' + 



2 10 nieirii 3 2 m? 5 



-» -I ^» 

3 M2 ^ 3M2'^ 



x/TOeTOi , 10 me „TOe „ 11 TOe - „ me ., 



, , 8 m^ 40 memi , Jmerrii , 16 me 16mi , , 4 m^ JrueXni 8 me 
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+ 



dcR^{De 



Jrrii 



1 me 



I mi 



P [Zei - \ Cp) + - — {l + Z{Zei - J C 

"me 2 M V me 



memt ^ jm[ 2^ , o "^^ y/mefnl jm[ juH 21 



'mi r„ ., I mi ,0 / mi . ^Jm^ml 



y/memi 



me 



I mi 



.rrii 



x(-6c(l - c2)p^^^^ - 2(1 - c'){l - ^c')^{zei - ^^cp) + (1 - 3c'){zei - a/^cp)[7 + 2^p 



+Acp 
+2A3 



x/memi , imi . , 



^m.em,i TOe 2n 



me 



M M 



pC (1 - C )(2;ei - 



cp) 



. , 7 me 7 mi o 21 m^. 



""'per) - A27py^{l - 3c'){Zei - .[^PC) 
mi M V TO,; 




) , 



524 



/■°° _ 2 r 47meTOi 5 2 5 to^ 2 2 TOemi 4 



6 Af2 3^2^ 



3M2' 



3 M2 



e-v^^ Imi m, 2 ,2me 2 4 to^ ^2x 



M 



3 M 



M 



3 M' 



5meM 



. , 8 TOeTOi 40meTOi 2 , 16 Jmemi mi , 2^^ , a mem,- Jmem^ 8 mi 

+A2P(--r^r7^ + T7^7^P +77 T7-^2(l+P )) + A3P(— - i? r TTTt) 



15 M2 15 M2 



15 M M 



15 M2 



dcRl{De 



y/memi y/memi 4 , „ Wg 3, 



M 



M 



-cp) 



+ 



m^mij2 
me ^ 



m, 

•-^'^p\l + Z{Zei - J^Cpf) + •^c'^PiZei ' J — Cp)[l + 3(^ei " J — Cpf] 

me M y me M V me V me 



me me 9 , , 
mi M 



M 15 M 

m. 



cp 



+A2 



[2 + 4(^ei - J—CPY + 3(^ei - W— CP)^]) 
M y TOe V me 

y/memi 



M 



[6c(-l + c^) v^p2 + 2(1 - c2)(l - 3c2)^p(ze. - ,/^cp) 



' mi 
me 



, (^ Q 2^ / x/o"^i 2 , 4 y/memi jnH o'^s/'i , ^ 

+ (1 - 3c )p{Zei - J — CP)(2^P + 4cp^^^(.e. - J — Cp) + 2^(1 + {Ze 



^m.em., m,,: 2/-, 2n . x 



) • 



The remaining expressions (i = 3, 4) follow as 

531 = bi3, (e ^ i)], 532 = [5i4,(e<-^i)], 533 = [511, (e i)]. 
541 = [523, (e ^ 542 = [524, (e ^ i)], 543 = [521 , (e i)] , 

For the collisions between identical species (e, i) we have 



534 = [512, (e i)], 
544 = [522, (e ^ i)]- 



1 eV 1 \ ' / 
n 



and 



A'^ = 1 



/rAv2 1 



2TOcfcT p2 

so that the contributions of electron-electron collisions = 1,2) follow as 

/ill = , 



hi2 = 4rJ j dp j dce-P'p{l - 3c^)[X''\n f ^T^j + 2](V2ze - cp)D{V2ze - cp) , 
h2i = Arl dp j dce-P' p{l - 3c2)[A" In (^^7^) + 2](\/2ze - cp)D{V2ze - cp) , 



(55) 



(56) 



poo f A 1 f\ 

h22 = dpe-P' [X^ In (j^) + 2] I 3/ + -r^p' + 2r^2 J ^ dcp 

X [(1 - 3c2)(%/2ze - cp)(p2 - p2^2 ^ g ^ 22^) + 3pc{c^ - l)]D{V2ze - cp)} . (57) 
The expressions for ion-ion collisions = 3,4) follow as 

h33 ^ [hii,{e ^ i)], = [^12, (e ^ i)]i /143 = [^21, (e ^ «)], /144 = [^22, (e ^ i)], (58) 

i.e. replacing the index e in cfj, A*^, Ze by the index i. 
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Figure Captions: 
Fig. 1: 

Real and imaginary part of the dielectric function as a function of the frequency at fixed wavenumber A: = 1 . 
The two-moment approach is compared with the one-moment approach and the RPA. 

Fig. 2: 

The same as Fig. 1 for wavenumber fc = 0.1 . 
Fig. 3: 

Imaginary part of the inverse dielectric function as a function of the frequency at fixed wavenumber k = 0.5 . 
The two-moment approach is compared with the RPA. 

Fig. 4: 

The same as Fig. 3 for wavenumber fc = 0.3 . 
Fig. 5: 

Imaginary part of the inverse dielectric function as a function of the frequency at different wavenumbers. 
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